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The random intensity of noise approach to one-dimensional Laval-DubruUe-Nazarenko type model 
having deductive support from the three-dimensional Navier-Stokes equation is used to describe La- 
grangian acceleration statistics of a fluid particle in developed turbulent flows. Intensity of additive 
noise and cross correlation between multiplicative and additive noises entering a nonlinear Langevin 
equation are assumed to depend on random velocity fluctuations in an exponential way. We use 
exact analytic result for the acceleration probability density function obtained as a stationary so- 
lution of the associated Fokker-Planck equation. We give a complete quantitative description of 
the available experimental data on conditional and unconditional acceleration statistics within the 
framework of a single model with a single set of fit parameters. The acceleration distribution and 
variance conditioned on Lagrangian velocity fluctuations, and the marginal distribution calculated 
by using independent Gaussian velocity statistics are found to be in a good agreement with the 
recent high-Reynolds-number Lagrangian experimental data. The fitted conditional mean accelera- 
tion is very small, that is in agreement with DNS, and increases for higher velocities but it departs 
from the experimental data, which exhibit anisotropy of the studied fiow. 
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I. INTRODUCTION 

Data analysis and modeling of statistical properties of a 
Lagrangian particle advected by a fully developed tur- 
bulent flow is of much practical interest and comple- 
ment traditional studies made in Eulerian framework. 
Strong and nonlocal character of Lagrangian particle 
coupling due to pressure effects makes the main ob- 
stacle to derive turbulence statistics from the Navier- 
Stokes equation. Recent breakthrough Lagrangian ex- 
periments Q, 1^ 1^ have motivated growing interest to 
a single-particle statistics. Some phenomenological ap- 
proaches 0, inspired by the non-extensive statistics 
formalism were used to describe Lagrangian accelera- 
tion of a fluid particle in a stationary developed turbulent 
flow within the framework of Langevin type equation; see 
also Refs. 7, 8, 9]. 

Some toy models of developed turbulence suffer from 
the lack of justification of a fit from turbulence dy- 
namics |14| . and the connection between specific non- 
thermodynamical processes and non-extensive mecha- 
nisms was argued to be generally not well defined [l5|. 
Recent one-dimensional (ID) stochastic particle mod- 
els and their refinements [lol ITlL were reviewed in 
Ref. |l3j | , in which importance of Navier-Stokes equation 
based approaches is emphasized. 

Fluid particle dynamics in a developed turbulent fiow 
is described in terms of a generalized Brownian motion 
with the Lagrangian acceleration of individual particle 
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viewed as a dynamical variable. In the data processing, 
the acceleration is associated with the Lagrangian veloc- 
ity increment in time for sufficiently small timescales, in 
a far dissipative subrange where turbulent fiuctuations 
are smoothed. Such models are generally based upon a 
hierarchy of characteristic timescales in the system and 
naturally employ one-point statistical description using 
Langevin type equation, or the associated Fokker-Planck 
equation for one-point probability density function of the 
variable. Noises entering the Langevin type equation are 
treated along a fiuid particle trajectory, and the Fokker- 
Planck approximation makes connection between the dy- 
namics and statistical approach. 

Recently we have shown jH El that the ID 

Laval-DubruUe-Nazarenko (LDN) toy model [T9ll20l| of 
the acceleration evolution with the model turbulent vis- 
cosity ft and coupled delta-correlated Gaussian multi- 
plicative and additive noises is in a good agreement with 
the high-precision Lagrangian experimental data on ac- 
celeration statistics [l|, |^ |^ ; Taylor microscale Reynolds 
number Rx = 690, the measured normalized acceleration 
range is —60 < a/(a^)^/^ < 60, and the resolution is 1/65 
of the Kolmogorov length scale of the fiow. The long- 
standing Heisenberg-Yaglom scaling of a component of 
Lagrangian acceleration, (a^) — aQV?l'^v~^l'^L~^^'^ , was 
confirmed experimentally yj] to a very high accuracy, for 
about seven orders of magnitude in the acceleration vari- 
ance, or two orders of the root mean square velocity m, at 
R\ > 500; flo is the Kolmogorov constant, i' is the kine- 
matic viscosity, and L is the Lagrangian integral length 
scale. Long-time correlations and the occurrence of very 
large fiuctuations at small scales dominate the motion of 
a fiuid particle, and this leads to a new dynamical picture 
of turbulence U^^. 
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The original 3D and ID LDN models were formu- 
lated both in the Lagrangian and Eulerian frameworks 
for small-scale velocity increments in time and space re- 
spectively. They are based on the Gabor transformation 
(Fourier transform in windows) and a stochastic kind 
of Batchelor-Proudman rapid distortion theory (RDT) 
approach to the incompressible 3D Navier-Stokes equa- 
tion 19], and thus have a deductive support from tur- 
bulence dynamics. A study based on direct numerical 
simulations (DNS) of the 3D LDN model in the regime 
of decaying turbulence has been made. 

The random intensity of noise (RIN) approach 
provides an extension of the ID LDN model viewed in the 
limit of small timescale r for which Lagrangian velocity 
increments are proportional to r: u{t + T) — u(t) = Ta{t). 

The main idea of the RIN approach is simply to ac- 
count for the recently established two well separated La- 
grangian autocorrelation time scales for the velocity in- 
crements 21] and assume that certain model parameters, 
such as intensity of noise, fluctuate at large timescale. 

An analysis of such a simple ID model can shed some 
light to properties of the 3D LDN model of Lagrangian 
dynamics. Recent development of the 3D LDN model 
can be found in Ref . , in which some new methods of 
turbulent dynamo problem have been exploited. 

The experimental data on the axial component of La- 
grangian acceleration a of polystyrene tracer particle in 
the Rx = 690 water flow generated between counter- 
rotating disks have been fitted by the stretched expo- 
nential il|,,2cii2jj 



P{a) = Cexp 



{l + \ha/b2r)bl 



(1) 



Here, h = 0.513 ± 0.003, 62 = 0.563 ± 0.02, and 
63 = 1.600 ± 0.003 are fit parameters, and C = 0.733 
is a normalization constant. The studied flow is highly 
anisotropic at large scales and this appears to affect small 
scales, resulting in a small skewness of the acceleration 
distribution and observable distinction in the distribu- 
tions of different components of the velocity. At large 
acceleration values, tails of the distribution ^ decay 
asymptotically as exp[— ]a]°''^] that implies a convergence 
of the fourth-order moment (a^) = /_oo o,'^P{0')da, as 
confirmed by the experiment with measured normalized 
acceleration values up to \a\ = 60. The fiatness fac- 
tor of the distribution which characterizes widen- 
ing of its tails (when compared to a Gaussian) is = 
{a^) / {a?)'^ = 55.1, that is in agreement the experimental 
value F = 55 ± 8. The Kolmogorov time of the fiow is 
— 0.93 ms. Low-pass filtering with the 0.23r^ width of 
the collected 1.7 x 10* data points was used, and the re- 
sponse time of the optically tracked 46 ^m tracer particle 
is 0.12t^. 

Recently, Chevillard et al. have constructed an 
appropriately recasted multifractal approach to describe 
statistics of Lagrangian velocity increments in a wide 
range of timescales, from the integral to dissipative one. 



The resulting theoretical distribution reproduces con- 
tinuous widening of the velocity increment probability 
density function (PDF) with the decrease of timescale, 
from a Gaussian-shaped to the stretched exponential 
as observed in Lagrangian experiments carried out at 
Cornell 0, S Q and ENS-Lyon [H |3, and DNS of 
the 3D Navier-Stokes equation. Two global parameters 
(Reynolds number and Lagrangian integral timescale) 
and two local parameters (intermittency parameter and 
smoothing parameter) with a parabolic singularity spec- 
trum were used to cover the data in the entire range of 
timescales. At dissipative timescale the obtained PDF 
fits the experimental data on Lagrangian acceleration to 
a good accuracy. The cumulant analysis made in this 
approach provides an understanding of the observed de- 
partures from the scaling when going from the integral 
to dissipative timescale. The used parabolic singular- 
ity spectrum D{h) is a hallmark of the log- normal (Kol- 
mogorov 1962) statistics and reproduces well the left- 
hand-side (corresponding to intense velocity increments) 
of the observed curve which is centered at 0.58 (> 1/2) 
but increasingly deviates at the right-hand-side of it (cor- 
responding to weak velocity increments). Another widely 
used statistics, the log-Poisson one, was shown to depart 
from Lagrangian observations in the same manner. The 
conditional acceleration statistics was not considered in 
this work. 

In a recent paper A. Reynolds |2^ developed a self- 
consistent second-order stochastic model with additive 
noise which accounts for dependence of the Lagrangian 
acceleration covariance matrix on Lagrangian velocities 
u. The observed dependence of the conditional acceler- 
ation variance {a?\u) on m was partially understood 
in terms of Lagrangian accelerations induced by vortex 
tubes within which the vorticity is constant and outside 
which the vorticity vanishes. Scaling relations were in- 
voked to derive a third-order polynomial structure of the 
isotropic covariance matrix as a function of squared ve- 
locity V? I23 . The inclusion of such conditional accelera- 
tion covariances in the model resulted in reduction of the 
predicted occurrence of small accelerations that meets 
experimental and DNS data for unconditional distribu- 
tions. The cores of the resulting conditional acceleration 
distributions were found to broaden with increasing u, in 
a qualitative agreement with the experiment. 

Sawford et al. |2^ have studied acceleration statistics 
from laboratory measurements and direct numerical sim- 
ulations in 3D turbulence at R\ ranging from 38 to 1000. 
For large \u\, the conditional acceleration covariance be- 
haves like . This is qualitatively consistent with the 
stretched exponential tails of the unconditional acceler- 
ation PDF. The conditional mean rate of change of the 
acceleration derived from the data has been shown con- 
sistent with the drift term in second-order Lagrangian 
stochastic models of turbulent transport. The correlation 
between the square of the acceleration and the square of 
the velocity has been shown small but not negligible. 

In very recent papers Biferale et al. [27I |2^ have pre- 
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sented interesting results of DNS of Lagrangian transport 
in homogeneous and isotropic turbulence with Rx up to 
280, a very accurate resolution of dissipative scales, and 
an integration time of about Lagrangian timescale. They 
have shown how the multifractal formalism offers an al- 
ternative approach which is rooted in the phenomenol- 
ogy of turbulence. The Lagrangian statistics was de- 
rived from the Eulerian statistics without introducing 
ad hoc hypotheses. Although the formalism is not ca- 
pable to account for small acceleration values (typical 
situation for the multifractal approach), the obtained ac- 
celeration PDF captures the DNS data well in the tails, 
with normalized acceleration values ranging from about 
|a|/(a2)i/2 = 1 up to |a|/(a2)i/2 = gg. Alas, one can ob- 
serve an overestimation in this range which can be clearly 
seen from the predicted contribution to the fourth-order 
moment, a^P(a), as compared to the DNS data. High 
degree of isotropy of the simulated stationary flow sug- 
gests equivalence of Cartesian components of accelera- 
tion aligned to fixed directions, and the resulting DNS 
distribution obtained by averaging over the components 
has been found with no observable asymmetry with re- 
spect to a ^ —a. The multifractal approach has been 
also used psj to obtain acceleration moments conditional 
on the velocity. Particularly, the multifractal prediction 
(a^|w) ~ u^-^^ agrees well with the DNS data for large 
velocity magnitudes. The predicted exponent 4.57 differs 
from the value 6 predicted recently by Sawford et al. 
and is very close to the Heisenberg-Yaglom scaling ex- 
ponent value 9/2. This indicates that the averaging of 
the above conditional acceleration variance {a^\u) over 
Gaussian distributed velocity u is consistent with the 
Heisenberg-Yaglom scaling law (see remark and Eq. (68) 
of Ref. Hi). 

In the present paper, we focus on ID LDN type dynam- 
ical modeling of the Lagrangian acceleration conditional 
on velocity fluctuations presented recently by Mordant, 
Crawford, and Bodenschatz in the experimental work Q . 
In contrast to our previous studies of the conditional ac- 
celeration statistics |0, 0, 0, 0| , here we give a self- 
consistent treatment of the model by explicit accounting 
for a Gaussian distribution of Lagrangian velocity fluc- 
tuations that is observed experimentally. We give a com- 
plete quantitative description of the available experimen- 
tal data on both the conditional and unconditional accel- 
eration statistics within the framework of a single model 
with a single set of fit parameters. Importance of the 
present approach is that the Lagrangian single-particle 
modeling is dynamical and has a deductive support from 
the Navier-Stokes equation, with few assumptions justi- 
fied by the turbulence phenomenology being used. This 
approach adds a new look to homogeneous isotropic tur- 
bulence modeling which is alternative to those given by 
the recent multifractal and ad hoc Langevin stochastic 
approaches. 

It should be emphasized that the Lagrangian veloc- 
ity is known to follow Gaussian distribution to a very 
good accuracy while the Lagrangian acceleration follows 



highly non-Gaussian distribution which is related to ex- 
tremely intermittent character of the acceleration, with 
pronounced central peak and relatively frequent acceler- 
ation bursts up to 80 standard deviations. We note that, 
theoretically, time derivative of a dynamical variable does 
not necessarily follow the same statistical distribution as 
that of the variable. 

Our consideration is restricted to a stationary one- 
point distribution function. Two-point statistical anal- 
ysis is of much interested and can be made elsewhere. 

The paper is organized as follows. In Sec. ^ we give 
a brief description of the ID LDN model and present 
the resulting acceleration distribution, which we treat as 
a conditional one by assigning stochastic properties to 
certain parameters. In Sec. IllII we briefly review results 
of our previous work and make sample fits of the obtained 
conditional and unconditional acceleration distributions 
and moments to the experimental data. In Sec. IIVI we 
discuss the obtained results and make conclusions. 



II. ID LAVAL-DUBRULLE-NAZARENKO 
MODEL OF SMALL-SCALE TURBULENCE 

In this Section, we present only a brief sketch of the 
ID LDN model and refer the reader to Refs. 0, for 
more details; see also Ref. This toy model can also 
be viewed as a passive scalar in a compressible ID flow. 

The main assumption of the LDN approach to the 
3D Navier-Stokes turbulence is to introduce and sepa- 
rate large-scale and small-scale parts in the 3D Navier- 
Stokes equation by using the Gabor transformation [T^ . 
This allows to consider analytically small-scale turbu- 
lence coupled to large-scale terms (the inter-scale cou- 
pling). The approach allows one to account for nonlocal 
interactions which were argued to be important in un- 
derstanding intermittency in developed turbulent flows. 
The other, large-scale, part of the equation can be treated 
separately (and, in principle, solved numerically given 
the forcing and boundary conditions) since the forcing is 
characterized by presumably narrow range of small wave 
numbers, and the small scales make little effect on it. 
Small-scale interactions are modeled by a turbulent vis- 
cosity and were shown numerically to make small contri- 
bution to the anomalous scaling (intermittency) in the 
decaying turbulence. Nevertheless, these are important 
when fitting model distribution to the experimental data. 
The 3D LDN model of small scale turbulence was used to 
formulate simplified ID LDN model, which was studied 
both in the Eulerian and Lagrangian frames . 

We use probability density function obtained as a sta- 
tionary solution of the Fokker-Planck equation that cor- 
responds to a consideration of statistically stationary 
state; statistical homogeneity and isotropy of the 3D flow 
is assumed as well. This equation is derived from the 
Langevin equation for a component of Lagrangian accel- 
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eration a{t) [11111, 



da 
'di 



(2) 



where Vt = -^/j/q + B'^a? /k? is the turbulent viscosity 
modehng small-scale interactions, i/q is constant kine- 
matic viscosity, B is free parameter measuring the con- 
tribution of nonlinearity in a to the turbulent viscosity, 
and k is wave number; dtk — —k^, k(0) — ko, to model 
the RDT stretching effect in one-dimensional case. 

In the original 3D LDN model based on the Navier- 
Stokes equation, ^{t) is related to the velocity derivative 
tensor and (T±{t) describes a forcing of small scales by 
large scales via the energy cascade mechanism (nonlocal 
inter-scale coupling). In the ID LDN model, these are 
approximated by a sufficiently simple statistics inspired 
by the Kraichnan ensemble used for turbulent passive 
scalar and the Kraichnan-Kazantsev model of turbulent 
dynamo: external Gaussian white-in-time noises along a 
fluid particle trajectory, 

m) - 0, mat')) - 2DS{t - t'), - 0, (3) 

(CTi(t)(7±(t')> = '2a5{t - t'), {i{t)a^{t')) = 2\5{t ~ t'). 

Here, I?, a, and A are free parameters measuring inten- 
sity of the noises and their cross-correlation, respectively. 
Bigger D and a means bigger contribution of the veloc- 
ity derivative tensor and the inter-scale coupling (both 
viewed here as short-time autocorrelated processes) to 
the small-scale dynamics. 

The model ^ puts an obvious limitation but is par- 
tially justified by DNS in the laboratory frame of ref- 
erence |l9j . The averaging is made over ensemble re- 
alizations. Zero means correspond to isotropy of the 
forces. Physically, the small scales are thus assumed to be 
stochastically distorted by much larger scales. We stress 
that a correlation between the noises ^ and a± is not ad 
hoc assumption but a consequence of their structure as 
they contain the same large-scale velocity serving as a 
unifying agent between the noises. 

It should be emphasized that the ID LDN toy model 
and its particular case Q have several limitations re- 
lated to the LDN separation of small and large scales 
allowing to study exclusively nonlocal effects associated 
to the linear process of distortions of small scales by a 
strain produced by large scales, the use of model turbu- 
lent viscosity, and one-dimensionality. 

In the Lagrangian frame the wave number k is replaced 
in terms of the initial value fcp = ^(0) and time while the 
parameters acquire dependence on fco [l9l |: we drop the 
subscript in fco in Eq. (0) and subsequent formulas to 
simplify notation. 

Thus, one makes a closure by treating the combined 
effect of large scales, for which one has a different dy- 
namical LDN equation that could be in principle solved 
numerically |23j |. and nonlocal inter-scale coupling, as a 
pair of given external noises. The large-scale dynamics 
is local in wave number space and hence it is weakly 



affected by the small scales. The price of the simplifica- 
tion (O is that one introduces free parameters to the de- 
scription. Matching small-scale dynamics to large-scale 
one deserves a separate study. Despite 3D turbulence 
is known to be more sensitive to large-scale forcing or 
boundary conditions, as compared to 2D one, the used 
simplification ijjjl is relevant for high-Reynolds-number 
flows to some extent 19., , 23], and allows one to advance in 
analytical treatment of the problem. It should be noted 
that scaling properties of the system described by Eq. Q 
reveal a robust character with respect to the selection of 
noises ^ and a± (see Ref. |0 and references therein). 

The acceleration PDF stemming from the stochastic 
model has been calculated exactly in our previous 

work El, 



P{a) = Cexp 



dx- 



-k'^x^vl+B'^x'^/B-Dx+X 



Dx"^ - 2\x + a 

Cexpi-i'tkyP + F{c) + F{~c)] 
{Da^-2Xa + ay/^{2Bka + ,y,k^)2B\k/D^ ' 



(4) 



for constant parameters. Here, C is normalization con- 
stant and we have denoted 



2D^ 



ci = 



2c2D^c ^ciC2ic-Da + \) 

X {B^{X^ + cA - Da)a + c(Di/f fc^ + C2i^t))], (5) 

c= -i^/Da- X^, (6) 

B2(4A3+4cA2-3£»aA - cDa)+D^{c+X)iy^k^ , (7) 



C2 



B2(2A2 + 2cA - Da)k'^ + D^v^k^. (8) 



The distribution Q is characterized by the presence of 
exponential cut off, complicated power-law dependence, 
and terms responsible for a skewness (asymmetry with 
respect to a ^ —a). 

One way when comparing the model with the exper- 
iment is to make a direct fit of the obtained PDF Q 
to the experimental data on unconditional acceleration 
distribution by assuming all the parameters and wave 
number to be constant. 

Particularly, this implies a reduction of the original ID 
LDN model since wave number is taken to be fixed so that 
the artificial ID compressibility aimed to model RDT 
stretching effect in ID case is not considered. We note 
that the Lagrangian acceleration is usually associated to 
the dissipative scale, and in the present paper we do not 
study dependence of the parameters on the wave number. 
Such a dependence for velocity increments was analyzed 
in Ref. 19] with the expected result that for larger scales 
the velocity increment PDF tends to a Gaussian form. 
The Gaussian form is reproduced also when D —^ and 
B — > 0, i.e., the process becomes purely additive with a 
linear drift term. 

Without loss of generality one can put, in a numerical 
study, fc = 1 and the additive noise intensity a = 1 by 
rescaling the multiplicative noise intensity Z) > 0, the 
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turbulent viscosity parameter B > 0, the kinematic vis- 
cosity uq > 0, and the cross correlation parameter A. 
The particular cases B = and i^o = at A = 0, and the 
general case at A = were studied in detail in Rcf. Ts'l. 
Nonzero A is responsible for an asymmetry of the PDF Q 
and in 3D picture corresponds to a correlation between 
stretching and vorticity (the energy cascade). Particu- 
larly, in the Eulerian framework the third-order moment 
of spatial velocity increment {{Siu)^) was found to be pro- 
portional to cross-correlation parameter, in accord to a 
kind of generalized Karman-Howarth relationship (l9l | . 

However, the approximation based on constant param- 
eters does not allow one to consider both the conditional 
and unconditional acceleration statistics. 

In the next Section, we extend the model ©-lEl by 
assuming certain model parameters in Eq. to be de- 
pendent on random velocity fluctuations. This extension 
is compatible with the 3D LDN approach as ^ and a± 
depend on velocity fluctuations and contain large-scale 
quantities due to their definitions [l^ . Such a functional 
dependence and longtime fluctuations have been ignored 
when making the simplification Q. We partially restore 
them. This is the main point of our consideration, and 
the functional form of the distribution is thus due to 
Eq. Q with certain parameters being now treated as 
functions of stochastic velocity u. Observations are that 
the acceleration variance does depend on the same com- 
ponent of velocity fluctuations. Local homogeneity as- 
sumed by Kolmogorov 1941 theory is thus broken that is 
a prerequisite to describe turbulence intermittency. The 
scaling approach indicates an essential character of such 
a dependence. Lagrangian intermittency is known to be 
much stronger than the Eulerian one due to existence of 
very intense vortical structures at small scales and ab- 
sence of the so called sweeping effect in the Lagrangian 
frame. 

We point out that characteristic time of variation of 
the parameters should be sufficiently large to justify ap- 
proximation that the resulting PDF @ is used with inde- 
pendent randomized parameters, P(a|Parameters). Two 
well separated timescales in the Lagrangian velocity in- 
crement autocorrelation have been established both by 
experiments and DNS . The large timescale has been 
found of the order of the Lagrangian integral scale and 
corresponds to a magnitude part that is in accord to our 
assumption that the intensity of noise along the trajec- 
tory is longtime fluctuating. 



III. THE CONDITIONAL ACCELERATION 
STATISTICS 

The experimental unconditional and conditional dis- 
tributions, which we denote for brevity by /'cxp(a) and 
Pcxpialu) respectively, were found to be approximately 
of the same stretched exponential form at u = (Fig.^l, 
and both reveal a strong Lagrangian turbulence intermit- 
tency (Jl- This similarity indicates that they share the 
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FIG. 1: A comparison of the experimental unconditional La- 
grangian acceleration PDF (dots) and the experimental con- 
ditional Lagrangian acceleration PDF at velocity fluctuations 
u — (line) PiQ; the acceleration component a is normalized 
to unit variance. 



same process underlying the intermittency. 
Accordingly, in our previous studies 

Ta, 111, [Hill 

we 

used the result of our direct fit of the PDF Q to Pexp(a), 
which was measured with a high precision; 3% relative 
uncertainty for \a\/{a^)^/^ < 10 2, 3]. We assumed that 
the parameters a and A entering Eq. depend on the 
amplitude of Lagrangian velocity fluctuations u, while 
D, B, and i^q are taken to be flxed at the fitted values 
{k — 1). Theoretically, only a and A depend explicitly 
on large-scale velocity due to 3D LDN model, while the 
other parameters not. 

An exponential form of a{u) has been proposed in 
Ref. and was found to be relevant from both the 
(Kolmogorov 1962) phenomenological and experimental 
points of view. Particularly, such a form leads to the 
log-normal RIN model when u is independent Gaussian 
distributed with zero mean j^, and yields the acceler- 
ation PDF whose low-probability tails are in agreement 
with experiments 6, 13] . Also, we used an exponential 
form of X{u) so that the conditional acceleration PDF 
Q takes the form P{a\u) = P{a\a{u), X(u)). Such a 
form was found to provide good fits of (i) the conditional 
probability density function P{a\u) to Pcxpia\u); (ii) the 
conditional acceleration variance (a^lu) ; and (iii) the con- 
ditional mean acceleration (oIm) |l7l | at various u that 
meet the experimental data. A brief report on these 
results is presented in Ref. [l8l |. 

However, a self-consistent consideration of the model 
assumes fitting of P(a\u) to Pcxp(a|i*), and the marginal 
PDF computed due to 

/oo 
P{a\u)g{u)du, (9) 
-oo 

where g{u) is PDF of independent velocity fiuctuations, 
should reproduce Poxp(a)- The marginal distribution cor- 
responds to a convolution of the stationary acceleration 
statistics with independent random velocity fiuctuations. 
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FIG. 2: Theoretical conditional acceleration PDF P{a\u) 
(line) and the experimental conditional acceleration PDF 
(dashed line) at velocity fluctuations u — 0, 0.45, 0.89, 1.3, 
1.8, 2.2, 2.7, 3.1 Q (from top to bottom, shifted by repeated 
factor 0.1 for clarity); the acceleration component a is nor- 
malized to unit variance, and the same component of velocity 
u is given in root mean square units. 
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FIG. 3: Theoretical conditional acceleration mean {a\u} (tri- 
angles) and the experimental conditional acceleration mean 
(squares) as functions of velocity fluctuations. 
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FIG. 4: Theoretical conditional acceleration variance {a^lu) 
(triangles) and the experimental conditional acceleration vari- 
ance (squares) as functions of velocity fluctuations. 



In the present paper, we fill this gap. Our task is 
to fit a variety of the experimental data, both on the 
conditional and unconditional statistics of acceleration, 
with a single set of fit parameters. For this purpose we 
use the following natural steps. 

First we fit P{a\u) = P{a\a{u), X{u)) given by Eq. Q 
to Pexp(fl|u) 13] assuming that the parameters depend on 
u in an exponential way, 

a{u) — aoexp[\u\/ua], X{u) — Xoe:sq)[\u\/u\]. (10) 

Hereafter, we use normalized acceleration a and velocity 
fluctuations u. The fit parameter set is D > 0, i^o > 0, 
B > 0, Ao, Ua > 0, and u\ > (ao — 1, k ~ 1). 
The relations in Eq. (|10() mean that the additive noise 
intensity and the correlation between the noises become 
higher for bigger velocity fluctuations \u\. 

We fit P{a\0) to Pcxp(a|0), that excludes Ua and u\ 
from consideration, by varying D, i^q, and _B at ao = 1 
and Ao = —0.005. We notice that the available condi- 
tional statistics Pexp(a|M) is low for high velocities, the 
presented acceleration range is small, — 14 < a < 14, 
so that a rather big uncertainty remains when determin- 



ing fit values of the parameters. Changes in shape of 
Pcxp(a|u) with u increasing from 7i = to u = 3.1 are 
captured independently by the fit parameters Ua and ux. 
The result is shown in Fig. |21 Good overlapping of each 
curve with data points at all fixed magnitudes of u has 
been achieved. 

Second we calculate the conditional mean {a\u) and 
the conditional variance (a^|u) and compare them with 
the experimental data. This decreases uncertainty in fit 
parameter values. The results are shown in Figs. 13 and 
21 Note that {a\u) as a function of u is very small that 
does not match the experiment. We will discuss this in 
Sec. Hvl below. 

Finally we calculate numerically the marginal distri- 
bution Pni{a) given by Eq. with the conditional PDF 
P{a\a{u), X(u)) and Gaussian distribution of velocity 
fluctuations, 



5(") = 



exp 



u 



(11) 



at fixed a ranging from —100 to 100 with the step 0.1. 
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FIG. 6: Lin-lin plot of the central part of the curves of Fig.|3 
Same notation as in Fig.|S] 



Then we make an interpolation and fit it to Poxp(a)- A 
noticeable effect of the integration over u with Gaussian 
g{u) is a widening of tails of the distribution that meets 
Fig. n the integration range — 20 < w < 20 has been 
used. The fit of Pm(a) to Poxp(a) strongly decreases the 
uncertainty but the most strict determination of fit values 
comes due to a comparison of the theoretical contribution 
to fourth-order moment, a*P{a), with the experimental 
data. The results are shown in Figs.|Sl|Sl and[7| Quahty 
of these sample fits is better than in the other recent 
stochastic models reviewed in Ref. [T^ . In particular, 
the core of the unconditional distribution reproduces very 
well that given by the stretched exponential as shown 
in Fig. 1^ However, both curves a bit underestimate the 
height at a = 0. 

The value Aq — —0.005 has been obtained by adjusting 
the theoretical curve to slightly different heights of the 
peaks of the observed a^P{a) shown in Fig.[7| Note that 
the model does not assume the use of ad hoc skewness of 
the forcing. Nonzero cross correlation parameter A natu- 
rally results not only in small mean acceleration but also 
in a skewness of both the theoretical distributions P{a\u) 



and Pm(a)- This skewness may be associated to the Eu- 
lerian downscale skewness generation, which despite of 
being small for homogeneous flows is known to be of a 
fundamental character in the inertial range (Kolmogorov 
four-fifths law), since the Eulerian {{Siu)^) was found to 
be proportional to cross-correlation parameter. 

We stress that the observed very small skewness of 
acceleration distribution is attributed to the effect of 
anisotropy of the studied flow. How the large-scale asym- 
metry affects smallest scales of the flow is an interesting 
problem. Our fit made by using nonzero A is of an il- 
lustrative character, to verify whether it can explain the 
observed increase of the conditional mean acceleration 
with increasing velocity depicted in Fig. |3| This issue 
will be discussed further in Sec. IIVI 

The following remarks are in order. Our finding is 
that the condition Ua < u\ provides a convergence of 
Pin(a)- Also, ux should not be small to provide assumed 
condition A <C a at arbitrary u (the cross correlation is 
small as compared to both noise intensities a and D) 0, 
IT^ . We used these criteria when making the fits. 

The resulting sample fit values are given by 

D = 2.1, j/o = 5.0, B ^ 0.35, (12) 
Ao = -0.005, Ua ^ 3.0, UA = 3.0, 

with ao = 1 and fc = 1. The theoretical curves in Figs.El- 
13 are shown for this sample set of values, which require 
a further fine tuning. Such a small value of A as com- 
pared to a or Z? is in agreement with that obtained in 
the LDN direct numerical simulations. The calculated 
flatness factor F = 49.3 of Pin(a) is in agreement with 
the experimental value 55 ± 8. 

To summarize, the considered Navier-Stokes equation 
based ID toy model l|H)- (|10|) is capable to fit all the 
available high-precision experimental data on the condi- 
tional and unconditional Lagrangian acceleration statis- 
tics Q, 1^ 13 with the single set of parameters H12|) to a 
good accuracy, with an exception being only the condi- 
tional mean acceleration. 
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IV. DISCUSSION AND CONCLUSIONS 

One can see from Fig. O that at the values of fit pa- 
rameters (|12f) the predicted conditional mean accelera- 
tion {a\u) qualitatively is in agreement but does not re- 
produce the experimental data. Namely, it is nonzero 
due to nonzero A and increases with the increase of |m| 
but remains to be very small even at high values of 
The conditional mean acceleration is evidently zero 
for a symmetrical distribution (A = 0) and should be 
zero for statistically homogeneous isotropic turbulence. 
The observed departure from zero is thought to reflect 
anisotropy of the studied flow albeit the DNS of homo- 
geneous isotropic turbulence also reveals slightly nonzero 
mean 0]. 

To reduce the discrepancy, we have tried the value 
u\ = 1.0 instead of u\ = 3.0 to provide faster increase of 
|A| for higher This implies a good fit to the experi- 
mental conditional mean acceleration (see, e.g.. Fig. 2 in 
Ref. 01) but we found an excess asymmetry of P{a\u) at 
high u, with big departure from observations, and diver- 
gencies when calculating (a) . The reason of the diver- 
gency is that X(u) at ux = 1.0 grows faster than a{u) at 
Ua = 3.0 so that A becomes comparable or bigger than a 
with increasing u, and when A^ — > Da the function F{c) 
defined by Eq. ||SJ) undergoes unbound growth. Thus we 
conclude that the observed conditional mean accelera- 
tion is mainly due to the flow anisotropy effect rather 
than some intrinsic dynamical mechanism associated to 
the developed turbulence. 

In general one observes a rather small relative increase 
of the conditional mean acceleration for higher \u\ that 
eventually reflects a coup ling of the acceleration to large 
scales of the studied flow |l9ll2^ . This coupling could be 
accounted also by introducing a correlation between the 
acceleration and velocity fluctuations. This possibility is 
of much interest to explore as it may yield the deficient 
increase of {a\u) but it is beyond the scope of the present 
formalism, which assumes an independent velocity statis- 
tics. We note also that in contrast to the experimental 
data on the variance (a^|u) the experimental (a|u) ex- 
hibits small asymmetry with respect to u — > — u (not 
shown in Fig.O)). 

In the present paper, the multiplicative noise intensity 
D was taken to be independent on the velocity fluctua- 
tions u. The effect of variation of D has been considered 
in Ref. 0| with the qualitative result that it does not 
provide the specific change in shape of P{a\u) observed 
in experiments. However, a weak dependence of Z? on u 
can not be ruled out. 

In summary, the presented ID LDN type stochastic toy 
model with the velocity-dependent additive noise inten- 
sity and cross correlation parameter is shown to capture 
main features of the observed conditional and uncondi- 
tional Lagrangian acceleration statistics to a good accu- 
racy except for the discrepancy in the conditional mean 
acceleration which can be attributed to certain coupling 
of the acceleration to large scales of the studied flow. 



The main result is of course not only good sample fits 
which are important to test performance of the model but 
also certain advance in understanding of the mechanism 
of Lagrangian intermittency provided by the dynamical 
Laval-DubruUe-Nazarenko approach to small-scale tur- 
bulence. 

The central point is that the LDN toy model has a 
strong deductive support from the Navier-Stokes turbu- 
lence. The obtained exact analytic result for the con- 
ditional acceleration distribution and the use of recent 
high-precision Lagrangian experimental data on condi- 
tional and unconditional acceleration statistics provide a 
detailed analysis of the mechanism within the adopted 
framework. Effects of large timescales (nonlocality) and 
turbulent viscosity (nonlinearity) have been found of 
much importance in Lagrangian acceleration steady-state 
statistics. The detailed study of conditional accelera- 
tion statistics have revealed a specific model structure 
of the external large-scale dynamics and nonlocal inter- 
scale coupling for homogeneous high-Reynolds-number 
flows. The additive noise associated to the downscale 
energy transfer mechanism encodes the main contribu- 
tion to the velocity dependence of the acceleration statis- 
tics. The cross correlation between the model additive 
and multiplicative noises associated to a correlation be- 
tween stretching and vorticity naturally provides a skew- 
ness of distributions and a nonzero mean. Weakness of 
this correlation is a theoretical requirement that meets 
the Lagrangian and Eulerian experiments and DNS of 
homogeneous isotropic turbulence. The observed con- 
ditional mean acceleration is mainly related to the flow 
anisotropy. The cross correlation is related to the four- 
fifths Kolmogorov law but the effect of skewness is neg- 
ligibly small as the result of relatively large intensity of 
the additive noise, which tends to symmetrize acceler- 
ation distributions. This is a dynamical evidence im- 
plied by the model rather than a direct consequence 
of a priori assumption on isotropy in the spirit of Kol- 
mogorov 1941 theory. The use of exponential depen- 
dence of certain noise parameters on statistically inde- 
pendent Gaussian distributed Lagrangian velocity fluc- 
tuations has been found appropriate to cover new experi- 
mental data on conditional statistics and to transfer from 
the conditional to unconditional acceleration distribution 
both exhibiting a strong Lagrangian intermittency of the 
flow. Such a dependence is also compatible with the log- 
normal statistics assumed by the Kolmogorov 1962 the- 
ory. The Gaussian white-in-time multiplicative noise and 
longtime correlated intensity of the additive noise were 
both found to make an essential contribution to intermit- 
tent bursts. 
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